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. ABSTRACT 

■ - - ' We present new measurements of the luminosity function (LF) of Luminous Red 

Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) and the 2dF-SDSS LRG 
and Quasar (2SLAQ) survey. We have carefully quantified, and corrected for, un- 
certainties in the K and evolutionary corrections, differences in the colour selection 
methods, and the effects of photometric errors, thus ensuring we are studying the same 
galaxy population in both surveys. Using a limited subset of 6326 SDSS LRGs (with 
0.17 < z < 0.24) and 1725 2SLAQ LRGs (with 0.5 < z < 0.6), for which the matching 
colour selection is most reliable, we find no evidence for any additional evolution in 
the LRG LF, over this redshift range, beyond that expected from a simple passive 
evolution model. This lack of additional evolution is quantified using the comoving lu- 
minosity density of SDSS and 2SLAQ LRGs, brighter than Mo. 2^ — 5logho,r = —22.5, 
which are 2.51 ± 0.03 x IO-^LqMpc"^ and 2.44 ± 0.15 x 10"'' Lq Mpc"^ respec- 
tively (< 10% uncertainty). We compare our LFs to the COMBO-17 data and find 
excellent agreement over the same redshift range. Together, these surveys show no 
evidence for additional evolution (beyond passive) in the LF of LRGs brighter than 
Mo.2r - 5logho,j = -21 (or brighter than ^ L*). We test our SDSS and 2SLAQ LFs 
against a simple "dry merger" model for the evolution of massive red galaxies and 
find that at least half of the LRGs at z ~ 0.2 must already have been well-assembled 
(with more than half their stellar mass) by z ~ 0.6. This limit is barely consistent 
with recent results from semi-analytical models of galaxy evolution. 

Key words: cosmology: observations - galaxies: abundances - galaxies: ellipticals 
and lenticular, cD - galaxies: evolution - galaxies: fundamental parameters 
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1 INTRODUCTION 

The formation of massive elliptical galaxies is a major 
conundrum in modern cosmology. Observationally, it is 
clear that such galaxies formed the bulk of their stars at 
redshifts greater than two, with evidence coming from a 
variety of studies, i ncluding observatio n s of clusters and 
groups of galaxies teower ct al.1 Il992t lE llis et alj Il997t 
Kodama et al. 1998; dc Propris et al.l Il999: Brough et al 



20021: IWakeet al] [20051: IHolden et alj|2005l IPimbblet et al 
200a), optical /NIR galaxy co unt and colour studies 
jMetcalfe et alJll996ll200ll^2005^ and spectr oscopic surveys 
of field galaxies over a range of redshifts fBernardi et alj 
1998: Tragcr ot al. 2000; Bernardi ct al. 2003 a b JdF 
Hogg et a l. 2002; Kuntschn ei^ fc DaviesI 119 98; Bal drv et alJ 
20041: iGlazebrook et aL 1120041: iMcCarthv et al.l 12004 



Cimatti et alJ l2004l: iThqmas et alJ 120051: I Papovich et alJ 
20051: iBernardi et alJl200^ 



These studies also suggest that the evolution of a 
majority of these massive ellipticals is consistent with a 
simple passive model of stellar evolution. For example, 
[Bernardi et al. ( 2003ft,3) have u sed the Sloan Digital Sky 
Survey fSDSS: IVork et al.ll200ol) to study the evolution of 
luminous early-type galaxies (out to z ~ 0.3) and find 
that, for a fixed velocity dispersion, older ellipticals are red- 
der and fainter in luminosity, fully consistent with the ex- 
pected fading of their stellar populations with time. They 
also find that the environment (as measured by local den- 
sity) of luminous ellipticals has only a weak effect upon 
their properties, e.g., the Fundamental Plane fades by only 
0.075 mags/arcsec^ from the field to the cores of clus- 
ters (consistent with cluster galaxies forming at a slightly 
earlier epoch than field ellipticals) . Several other studies 
have also found very similar results ijK imtsch ner et alJl2002l : 
IClemens et~all l2006l iBernardi et al.ll2006l) . However, sev- 
eral authors have found evidence for recent and/or on- 
going star-formation in a fraction of l o cal massive ellipticals 
jTrager et alJl2000l : iGoto et aljl2003l : iFukugita et"aLl 12004 
iBalogh et alJl2005^. This fraction appears to increase with 
redshift" ilL^Borgne et_ad] j^005: Roscboom ct al. 2006) and 
decrease with mass JCaldwell et alJl2003l : iNelan et al.ll2005t 
IClemens ^8111200^) . 

One potential problem for studies of massive 
elliptical/early- type galaxies, particularly those track- 
ing the evolution with redshift, is that of progenitor bias 
Jvan Dokkum fc Franxlll996l: Ivan Dokkum et alJl2000l) . By 
studying only galaxies that look like massive early-type 
systems at high redshift, some of the progenitors of the 
low redshift massive early-type galaxies may be missed. 
An effective way to counter this problem is to study the 
evolution in the number density, or luminosity function 
(LF), of massive early-type galaxies with redshift. If these 
galaxies did indeed form at high redshift, then there will be 
little evidence for the evolution of their LF with redshift. 
Such studies have only just begun, e.g., the COMBO-17 
(C17) survey, which reported on the evolution of ~ 5000 red 
(and thus implied early-type) galaxies and found at most 
a factor of two evolutio n out to z 1 ( corresponding to 9 
Gyrs in look-back time: lBell et al.ll2004^ . This is consistent 
with earlier (but smaller in number) studi es of high redshift 
red (early-type) ga laxies from t he CFRS (Lill v et alJll995l : 
ISchade et akl Il999^ . CNOC2 (iLin et al.l 119991) and K20 



JPozzetti et alJl2003^ surveys. At higher redshifts {z > 1), 
recent spectroscopic surveys suggest that a significant 
fraction of massive galaxies arc alrea dy in place at the se 
early epochs ( jGlazebrook et al. 2004 : Cimatti et al.ll2004^ . 

Theoretically, the existence of massive, passively- 
evolved ellipticals has been a major challenge for some mod- 
els of galaxy evolution. For example, in the favoured hier- 
archical Cold Dark Matter (CDM) model of structure for- 
mation, such massive galaxies are expected to reside in 
the largest dark matter haloes, that form at late times 
through the merger of smaller mass haloes (see iKauffmanrJ 
(jl996) for early predictions). In recent years, there have 
been a number of prescriptions proposed to solve the "anti- 
hierarchical" nature of the formation and evolution of mas- 
sive ellipticals and thus better match the observations dis- 
cussed above, including: i) The reduction of the gas cooling 
rate (^Benson ct al. 2003); ii) Super- winds tha t eject the gas 
once i t has cooled, but be fore it can form stars iBenson et alJ 
I2OO3I ; ISaugh et 811120051) ; Hi) Shock heating of infalling gas 
and -PdV work of the gas iNaab e^ lj|200 ^: iv) Feedback 
from Active Galactic Nuclei fAGN: lKawata fc Gibsoiill2005l : 
[Scannapieco et al. 2005). 

The latter of these proposed effects (AGN feedback) 
is appealing because of the observed empirical relationship 
between the luminosity, central concentration and mass of 
galactic bul ges and the m ass of their central super-massive 
black holes iFerrarese fc Merritt 2000,;i'&emaine et al.l2002l: 



iNovak et al.ll2006l) . AGN feedback comes in two flavors 
i) The merger of gas-rich galaxies causing an initial star- 
burst, followed by the growth of the central black hole 
and a "quasar wind" which quenches further star-formation 
iHopkins et alj l2005l) ; ii ) Radio feedback from low lumi- 
nosity AGN that suppresses the cooling of gas in mas- 
sive halos resulting in the termination of late-time star- 
formation (Bower ct al. 2005; Croton ct al. 2006). Detailed 
simulations of these AGN feedback models on the evolu- 
tion of ellipticals has helped resolve the discrepancy between 
present observations a nd the naive hierarch i cal expectations 
of the CDM model fSpringel et al.' '2005'; 'De Lucia ct ajJ 
2005; Bower ct al. 2005; Hopkins etal,.,,20 05; Croton et al] 
[20061) . For example. [De Lucia etalJ I2OO5I) predict that for 
galaxies more massive than 10^^ M©, over 50% of their stars 
have formed by a median redshift oi z = 2.6, while the typ- 
ical assembly redshift of these galaxies (when these stars 
reside in a single object) is only z ~ 0.8. This requires "dry 
mergers" of the smaller galaxies (i.e. without gas) to avoid 
causing bursts of new star-formation. S ever al authors have 
found evidence for such dry mergers f see [van Dokkum 20o|; 
HeireraOliooi). 

In this paper, we present an accurate measurement 
of the evolution of the lu minosity function (LF) of Lumi- 
nous Red Galaxies (LRGs: fEisenstein et al.[[200l[) . that are 
predominantly massive elliptical galaxies and thus provide 
strong constraints on the physical mechanisms that have 
been suggested to produce anti-hierarchical galaxy forma- 
tion in a CDM universe. We utilise two samples of LRGs, 
one from the SDSS (0.15 < z < 0.37), and a new survey of 
higher redshift LRGs (0.45 < z < 0.8) selected from the 
multi-colour SDSS photometry, but spectroscopically ob- 
served using the 2dF spectrograph on the Anglo-Australian 
Telescope (A AT). This new survey is known as the 2dF- 
SDSS LRG and QSO (2SLAQ) survey. 
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The LRG component of the 2SLAQ survey is a signif- 
icant advance over previous surveys of massive early-type 
galaxies due to an increase in both the volume surveyed 
(~ 10^ Mpc''), thus reducing the problem of cosmic vari- 
ance, and the number of galaxy redshifts obtained, thus 
greatly increasing the statistical accuracy of the LF at bright 
magnitudes. In Section |5| we describe the 2SLAQ survey, 
while in Section 3 we present a detailed discussion of the K 
and evolutionary corrections, photometric errors and consis- 
tent colour selections for the SDSS and 2SLAQ samples. In 
Section r4.2l we present the luminosity functions of the SDSS 
and 2SLAQ surveys. In Section |3 we discuss our results in 
terms of models of galaxy evolution and conclude. Through- 
out, we adopt a cosmology of = 0.3, = 0.7 and Hq 
= 70 km s"^ Mpc"^ 



pleteness and uniform selection. The overall success rate of 
obtaining redshifts from the 2dF spectra for Sample 8 LRGs 
is 95%, while the centers of the 2dF fields were spaced by 
1.2°, resulting in an overall redshift completeness of sam- 
ple 8 LRG targets of ~75% across the whole survey area 
(see section 14. H . These data have re cently been used to 
calib r ate LRG photometri c redshifts llPadmanabhan et alJ 
120051: ICollister et al. llioo^) . 

Although the SDSS magnitude system wa s designed 
to be on the AB scale iFukugita et alJ Il996tl . the final 
calibration has differences from the proposed values by 
a few percent. We have applied the corrections viab = 
rasDSS + [-0.036,0.012,0.010,0.028,0.040] for u, g, r, i, z 
respectively (Eisenstein, priv. comm.). All magnitudes and 
colours presented t hroughout this pape r are corrected for 
Galactic extinction JSchlegel et alJll998h . 



2 LRGS IN THE SDSS AND 2SLAQ SURVEYS 

We present in this paper an analysis of LRGs taken from 
both the SDSS survey and the 2dF-SDSS LRG and QSO 
(2SLAQ) survey. For the SDSS data, we only use the LRG 
Cut I photometric selection (using the g — r and r — i 
colours of galaxies) as defined and discussed in Eisenstein 
et al. (2001; EOl). This provides us with a pseudo volume- 
limited sample of LRGs, with Mr ^ —21.8 and 0.15 < 
z < 0.35, selected from the SDSS Data Release (DR3, 
lAdelman-McCarthv et alJ l2005l) dataset. Below z = 0.15, 
the space density of SDSS LRGs increases by nearly 50% be- 
cause of contamination by low redshift star-forming galax- 
ies in the colour selection, while above z > 0.35, the space 
density of SDSS LRGs begins to decrease because of the fiux 
limit and the degeneracy between the colours and redshift of 
LRGs close to z ~ 0.37 where the 4000 Abreak feature passes 
out of the SDSS g-band into the ) — band ||Fukueita ct alj 
[1996). Within the range 0.15 < z < 0.35, the space density 
of SDSS Cut I LRGs is approximately constant with redshift 
(see EOl). 

In 2003, the 2SLAQ survey began with the goal of pro- 
ducing a pseudo volume-limited sample of 10,000 LRGs, 
with a median redshift of z = 0.55, and 10,000 faint z < 3 
quasars, both selected from the SDSS multi-colour imaging 
data. In this paper, we focus on the 2SLAQ LRGs which are 
selected using similar criteria as the Cut II SDSS LRGs in 
EOl. The key difi'erences are: (i) The apparent magnitude 
limit has been lowered to (model) < 19.8, thus extend- 
ing the volume-limited LRG samples to z ~ 0.6; (ii) The 
effective rest-frame colour cut {c± in EOl) has been shifted 
slightly bluer than in the SDSS LRG selection to accom- 
modate the density of 2dF fibres available. This provides a 
less conservative colour cut, at these higher redshifts, which 
is essential for studying potentially small changes in their 
colour. 

The details of the 2SLAQ LRG selection and o bserva- 
tions are presented elsewhere ijCannon et al.ll2006^ . How- 
ever, we have measured over 11000 LRG redshifts, covering 
180deg^ of SDSS imaging data, from 87 allocated nights of 
AAT time. Over 90% of these galaxies are within the range 
0.45 < z < 0.7. The tar geted LRGs where s plit into three 
subsamples as detailed in lCannon et alj ll2006f) . with the pri- 
mary sample (Sample 8) accounting for two thirds of these. 
We only focus on Sample 8 in this paper due to its high com- 



3 SAMPLE SELECTION 

The galaxy samples discussed in this paper were originally 
selected under the assumption that LRGs are old, passively- 
evolving galaxies. Here, we continue this fundamental as- 
sumption when investigating the evolution of the LRG lu- 
minosity function by applying the same passively-evolving 
models to our data when computing K and evolutionary 
(e) corrections, luminosity functions and colour selection 
boundaries. If this assumption is correct, then our K-l-e cor- 
rected luminosity functions will be identical over the whole 
redshift range studied here, with any additional evolution in 
the LRG population (beyond the simple passive model) be- 
ing displayed as a change in the LF with redshift. However, 
as we are using two slightly different selections of LRGs to 
generate our total LRG sample, it is vital that we select the 
same galaxy population at all redshifts. 

3.1 K+e corrections 

In order to make a fair comparison between the SDSS and 
2SLAQ LRG samples (at different redshifts), we must cor- 
rect the properties of our observed galaxies (magnitude, 
colours etc.) into rest-frame quantities by applying K- 
corrections. In addition, we can also correct these rest-frame 
quantities for the expected evolutionary changes over the 
redshift range studied. These e-corrections are usually per- 
formed assuming a model for the galaxy spectral energy dis- 
tribution (SED) and its evolution with redshift. In this pa- 
per, we therefore generat e our K-l-e corrections using the 
iBruzual fc CharlotI ^2003^ stellar population synthesis code. 
In detail, we generate two stellar population models the first 
of which forms all its stars in a single instantaneous burst at 
z = 9.84 (solar metallicity) and then evolves passively with 
no further star formation. The second model forms the bulk 
of its stars in a similar burst, but includes a small amount of 
continuous star formation throughout the rest of its evolu- 
tion, accounting for 5% of its final mass. These two models 
are shown in Figure 1 and are labelled "Passive" and "Pas- 
sive-|-SF" respectively. 

Figure0shows the colour evolution as a function of red- 
shift for these two models along with the actual measured 
colours of the SDSS and 2SLAQ LRGs. Some differences be- 
tween the models and data are evident, e.g., the models are 
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Figure 1. The g - r and r - i colours of the SDSS and 2SLAQ 
LRGs as a function of redshift (points) . The solid black line shows 
the median of the fits with the coloured lines showing the model 
tracks described in the text with (solid) and without (dashed) 
evolution applied. 



too red in g — r and too blue in r — i for the lowest redshift 
LRGs, with the opposite effect for the highest redshift LRGs. 
This offset between models and data in the lower redshift 
LRGs was noted in EOl and a correction to the g — r colours 
of 0.08 magnitudes was applied. However, with the addition 
of the higher redshift sample, it is clear that a simple offset 
is an inadequate correction over our entire redshift range. 
Such differences between the models and obse rved colours 
of earl y-type galaxies has been seen before, e.g.. lWake et alJ 
||200^ were unable to match the red sequence colours in 
galaxy clusters at similar redshifts, while a similar offset 
is seen in the colour-redshift plots of iFerreras et al.) ll2005l) 
from the GOODS/CDFS fields. Simple changes to the mod- 
els, such as the formation redshift, metallicity and IMF of 
the models, are unable to improve the model fits to the data. 
We also note t hat using the PEGASE stellar p opulation 
synthesis model jFioc fc Rocca-Volmerang^ll^97^ produced 
very similar results. It appears that the models are not accu- 
rately reproducing the shape of the spectrum between 4000A 
and 5000A causing an offset in t; — r and r — i as the r fil- 
ter passes through this (rest-frame) wavelength region. The 
g — i colours of LRGs are well reproduced by the models. 

In order to minimise the systematic uncertainties in the 
models, and thus uncertainties in our K-l-e corrections, we 
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Figure 2. The spectral energy distribution of the passive LRG 
model described in the text, along with the filter transmission 
curves blueshifted to 2 = 0.2 and z = 0.55. 



will exploit the fact that the redshifted SDSS g, r, i and 
z passbands ai z — 0.55 (i.e., the median redshift of the 
2SLAQ LRGs) approximately overlap the u, g, r and i pass- 
bands at z = 0.2 (the typical redshift of a SDSS LRG). This 
effect is illustrated in Figure |5| where the transmission curves 
of the SDSS filters are plotted a,t z = 0.2 and z — 0.55 and 
compared to the rest -frame spectral energy distribution of 
our passively evolved stellar population model as discussed 
above. 

Therefore, for each LRG, we estimate its K-fe correc- 
tion by interpolating between the two SED models using the 
observed g — i colour of the galaxy. For the SDSS LRG sam- 
ple, we simply correct the observed g, r and i magnitudes to 
g, r and i a.t z = 0.2. The median size of these K-|-e correc- 
tions ranges from 0.01 in i passband to 0.43 in g-band. For 
the 2SLAQ LRG sample, we again use the observed g — i 
colour to interpolate between the two SED models, but this 
time correct from observed r, i and z passbands to g, r and i 
at z = 0.2. The median size of these K and K-l-e corrections 
are 0.32 and 0.7 respectively when going from the r-band 
to the g-band, 0.1 and 0.4 going i-band to r-band, and 0.1 
and 0.4 when correcting from 2-band to i-band. Tables of 
these corrections for the two different SED models are given 
in Appendix 1X1 

Throughout the rest of this paper, we use the notation 
Mo. 2^ to represent the absolute magnitude of an LRG ob- 
served through a SDSS r-band fil ter redshifted to 0.2. This 
notation is similar to that used in lBlanton et alJ i2003l) . ex- 
cept here we include an evolutionary correction in addition 
to the K-correction. We note that Mo. 2^ ~ Mr(z = 0) -|- 0.11 
for the colour of a typical LRG. 



3.2 Consistent Colour Selection 

The differences in the SDSS and 2SLAQ selection cri- 
teria discussed above are explicitly illustrated in Figure 
El Here, we have t aken a sample of SDSS main galaxies 
Istrauss et alj|2002l) with 0.1 < 2; < 0.2 and K-l-e corrected 
all of them to 2 = 0.2. The SDSS main sample consists of all 
galaxies with r < 17.77 and provides a representative sample 
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Figure 3. The ^-^(g - i) versus Mo. 2^ colour magnitude relation 
for SDSS main galaxies with 0.1 < 2 < 0.2 all K+e corrected 
to z = 0.2. The small black points in each panel show the whole 
sample. The second panel shows those galaxies that would have 
been selected by the 2SLAQ selection criteria when K+e corrected 
to z = 0.55 (large points). The third panel shows those galaxies 
that would be selected by the SDSS LRG Cut I selection criteria 
when K+e corrected to 2 = 0.2 (large points). The final panel 
shows the galaxies that satisfy both the SDSS and 2SLAQ LRG 
selection criteria at both 2 = 0.2 and 0.55 respectively (large 
points). 



of the whole galaxy population. For this sample, we gener- 
ate the K+e corrections by interpolating between a passive 
model and one with continuous star formation based on the 
observed g ~ i colour of the galaxy. We then applied the 
SDSS and 2SLAQ LRG selection criteria at the two redshifts 
respectively. We plot in Figure |21 the colour-magnitude re- 
lation {^'^{g ~ i) versus Mo.2r) for these galaxies, where ^'^g 



is the g filter at z = 0.2, illustrating which galaxies would 
be selected by each criteria and their combination. This fig- 
ure clearly illustrates the bluer 2SLAQ selection as well as 
the magnitude dependence of the SDSS LRG colour selec- 
tion (EOl) as the selection cuts through the red sequence 
starting at Mo. 2^ > -22.8. 

Although the detailed colour evolution of LRGs remains 
unknown, we will proceed by assuming a simple passive 
model for LRGs evolution and thus make a self-consistent 
colour selection at both z — 0.2 and z = 0.55. We can then 
check whether the observed LF evolution is consistent with 
this simple hypothesis and whether any further evolution 
(beyond passive) is required to explain our observations. To 
do this, we first use the K+e corrections described above to 
correct all the LRGs in both samples to the redshift of the 
other sample i.e. we correct the 2SLAQ LRGs to z — 0.2 
and the SDSS LRGs to 2: = 0.55. We then require that for 
any individual LRG to be included in our analysis of the 
LRG LF it must satisfy both the SDSS criteria (at z = 0.2) 
and the 2SLAQ criteria (at z = 0.55). As might be expected, 
considering the broader colour and magnitude ranges of the 
2SLAQ selection criteria, most of the SDSS LRGs (~90%) 
would still be selected as LRGs a,t z — 0.55. The opposite 
is not true with only ~30% of the 2SLAQ LRGs satisfying 
the stricter SDSS Cut I LRG criteria. 

Figure|3]shows the colour distributions of all the 2SLAQ 
LRGs (red) and all the SDSS LRGs (blue) convolved with 
the typical photometric errors of the 2SLAQ LRGs, K+e 
corrected to z = 0.2. The left-hand panel shows the colour 
distributions of the full sample of SDSS and 2SLAQ LRGs 
considered in this paper (i.e., the raw data), while the middle 
panel shows only those galaxies which satisfy both the SDSS 
and 2SLAQ LRG selection criteria at z = 0.2 and 2 — 0.55 re- 
spectively. As mentioned previously, the 2SLAQ LRGs in the 
left-hand panel have bluer colour distributions compared to 
the SDSS LRGs, reflecting their bluer and fainter selection 
criteria. When the joint SDSS and 2SLAQ selection criteria 
are applied, the colour distributions of this restricted set of 
2SLAQ LRGs now become redder and much closer to the 
SDSS LRG colour distributions. However, an offset between 
the 2SLAQ and SDSS LRG colour distributions is still evi- 
dent (in the middle panel), i.e., 0.05 magnitudes in the mean 
g ~ r colour of LRGs. We believe this offset is due to resid- 
ual inaccuracies in the models used for our K-corrections, 
although we have tried to reduce this problem by correcting 
between overlapping filters. The wide redshift range of our 
samples (particularly in the SDSS) results in the errors in 
our K-corrections still being significant and detectable. 

To improve the agreement between the 2SLAQ and 
SDSS LRG samples, we restrict the redshift ranges of these 
two samples to be closer to the redshifts where the SDSS 
filters have their greatest overlap (Figure 2). We restrict 
the SDSS LRG sample to 0.17 < 2 < 0.24 and 2SLAQ 
LRG sample to 0.5 < z < 0.6, which reduces the ampli- 
tude of the K-corrections a.t z — 0.2 for the SDSS LRGs 
to a range of 0.005 to 0.06, compared to 0.01 to 0.11 with- 
out this restricted redshift range. Likewise, the amplitude of 
the 2SLAQ LRG K-corrections reduce to a range of 0.08 to 
0.30, compared to 0.11 to 0.36 without the restricted red- 
shift range. The right-hand panel of Figure 3] shows the 
colour distributions of LRGs that satisfy both the SDSS 
and 2SLAQ selection criteria within the restricted redshift 
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Table 1. The median photometric errors for the SDSS and 
2SLAQ LRGs from the single— epoch SDSS photometry and the 
multi-epoch SDSS photometry described in the text. 



u 


9 


r 


I 


z 


SDSS single-epoch 0.56 


0.04 


0.01 


0.01 




2SLAQ single-epoch 


0.15 


0.05 


0.03 


0.09 


2SLAQ multi-epoch 


0.05 


0.02 


0.01 


0.03 



ranges discussed above. This results in the colour distribu- 
tions of the two LRG samples being almost identical; the 
median colour difference is now less than 0.01 magnitudes. 
These restricted redshift LRG samples give us greater con- 
fidence that, under the assumption of passive evolution, we 
are now selecting the same type of galaxy in both the SDSS 
and 2SLAQ surveys. It also suggests that the discrepancy 
between the observed and model colours is only affecting the 
K-corrections and not the evolution corrections. The final 
redshift -restricted samples contain 6326 SDSS LRGs, with 
0.17 < z < 0.24, and 1725 2SLAQ LRGs, with 0.5 < z < 0.6. 
Clearly making such tight redshift cuts has resulted in a sig- 
nificant reduction in the size of our samples. However, we 
believe that minimising the errors in the K-l-e corrections is 
vital to ensure that the LRG samples from the SDSS and 
2SLAQ are as close as possible. 

3.3 Photometric errors 

Another potential bias that could affect the sample selec- 
tion is the larger photometric errors on the 2SLAQ LRGs 
compared to the SDSS LRGs (because they are intrinsically 
fainter and from the same imaging dataset). This could sys- 
tematically change the selection, as a function of redshift, 
as more 2SLAQ LRGs could be scattered both in and out 
of the sample as the photometric error increases. 

We attempt to measure this effect by utilising the multi- 
epoch SDSS imaging data available over a subsamp le of the 
2SLAQ survey area (see .Baldrv et al, 2005: Scranton et alJ 
I2OO5I) . This multi-epoch data covers a total of ~190 deg^ 
of the southern part of the 2SLAQ survey, and provides 
better signal-to-noise photometry for LRGs in this area as 
demonstrated by Tabled We begin by comparing the num- 
ber of LRGs that satisfy the 2SLAQ selection criteria using 
both the single and multi-epoch photometry. We find that 
10059 2SLAQ LRG targets were selected from the single- 
epoch photometry compared to 10265 2SLAQ LRG targets 
selected using the multi-epoch data. However, 25% of these 
targets are different between the two samples and have been 
scattered across the selection boundaries in almost equal 
numbers. In Figure |H| we show the colour-magnitude rela- 
tionship for LRGs both in common between the multi and 
single-epoch photometry as well as the LRGs only selected in 
one of the two data sets. As expected, LRGs selected only in 
the single-epoch data (i.e., missed by the 2SLAQ selection 
using the multi-epoch photometry) are fainter and bluer 
than 2SLAQ LRGs scattered out of the single-epoch selec- 
tion but selected using the multi-epoch photometry (the 
right-hand panel of Figure 5). We note however that the 
colour-magnitude relationship of LRGs for both the multi 
and single-epoch photometry is nearly identical for i < 19.3. 



This magnitude corresponds to Mo. 2^ = -23.0 at z = 0.6, the 
upper redshift limit used herein to calculate the luminosity 
function. 

As a further test, we can limit this comparison, to only 
LRGs with measured redshifts. This is possible because a 
subset of the 2dF fibres (~ 30%) were allocated to galaxy 
targets that lie slightly beyond the original 2SLAQ colour 
selection boundaries for the highest pri ority "Sample 8" se- 
lection (for details seelCannon et al. 2006). These extra LRG 
redshifts allow us to investigate the effect of photometric er- 
rors on the completeness of "Sample 8", although we are 
limited by the small size of the colour region beyond sample 
8 when considering those galaxies that weren't selected but 
should have been. In Figure|S|(top) , we show the multi-epoch 
colour and magnitude distributions for 2SLAQ LRGs, in the 
redshift range 0.5 < z < 0.6, all K-l-e corrected to z = 0.2 
as before. As above, 22% of the 2SLAQ LRGs are scattered 
into the sample because of photometric errors, i.e., they sat- 
isfy the selection criteria in the single-epoch photometry, 
but fail the criteria for the multi-epoch photometry. Again, 
these LRGs are bluer and fainter (in absolute magnitude) 
as shown by the dotted blue line in Figure |S| (top). How- 
ever, brighter than Mo. 2^ < = -22.65 the magnitude dis- 
tributions of those selected by the single and multi-epoch 
data are almost identical. We are therefore confident that 
the photometric errors have minimal effect on the 2SLAQ 
sample function brighter than this limit. 

A more significant effect of the photometric errors oc- 
curs when we apply the additional selection criteria we use 
to match the samples as discussed in Section 13.21 We can 
again try to quantify this effect using the multi-epoch data 
and again show in Figure El (centre) the magnitude distribu- 
tions of the single and multi-epoch selected 2SLAQ LRGs. 
However, here we only show those LRGs in the redshift range 
0.5 < z < 0.6 that pass both both the SDSS and 2SLAQ se- 
lection criteria. Unlike previously, where we had a limited 
set of spectra for galaxies beyond the selection boundaries, 
in this instance we have many galaxies far beyond the selec- 
tion boundaries due to the much bluer and fainter 2SLAQ 
selection. This is clearly visible when comparing the second 
and fourth panels of Figure |21 On inspection of the magni- 
tude distributions in Figure |S] it is clear that a significant 
fraction of galaxies are being scattered across the selection 
boundaries even for LRGs as bright as Mo. 2,. — -23.3. This 
is more clearly illustrated by the bottom panel of Figure 
1^ where we plot the ratio of the number of LRGs selected 
using the multi-epoch data to the number selected using 
the single-epoch data as a function of absolute magnitude. 
Since we are confident that we have fully sampled the colour- 
magnitude space beyond the selection boundaries we can use 
this ratio to correct the luminosity function and integrated 
number and luminosity densities presented in the Section 
14.21 We note that the only place that we are not sampling 
beyond the boundary is for the faintest and reddest objects 
and we are thus not confident of the correction and any re- 
sulting quantities for Mo.2r > -22.4. In order to generate 
accurate errors on this correction, we take the multi-epoch 
selected sample and add random Gaussian errors typical of 
the 2SLAQ photometric errors to each magnitude. We then 
calculate how many LRGs would be selected in our final 
sample. We repeat this procedure 10000 times and measure 
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Figure 4. The g — r and r — i colours of all the SDSS and 2SLAQ LRGs K+e corrected to 2: = 0.2. The SDSS distributions have been 
convolved with the typical 2SLAQ photometric errors. The left-hand plot shows all the LRGs in these two samples, while the middle 
plot shows only those LRGs (in both samples) that match both the SDSS and 2SLAQ LRG selection criteria (described in the text). 
The right— hand plot shows those LRGs that match both the SDSS and 2SLAQ LRG selection criteria and have an additional redshift 
restriction; 0.17 < z < 0.24 for the SDSS LRGs and 0.5 < z < 0.6 for the 2SLAQ LRGs. 



the standard deviations in each bin which are shown as the 
errors in Figure |^ (bottom). 

So far we have only discussed the effect of the photo- 
metric errors on the 2SLAQ LRGs. It is also worth briefly 
discussing any potential effect on the SDSS LRGs. As listed 
in Table0the typical photometric errors on the SDSS LRGs 
are much smaller than those of 2SLAQ and comparable to 
the multi-epoch errors on the 2SLAQ LRGs, except in the 
case of the u band. However, the u band data are only used 
when applying the 2SLAQ selection criteria to the SDSS 
LRGs and since the 2SLAQ criteria are typically signifi- 
cantly bluer and fainter than the SDSS criteria one would 
expect a small effect. In fact only 5% of the SDSS LRGs 
are removed when the 2SLAQ criteria are applied to them 
and most of these (3%) are as a result of the 2SLAQ mag- 
nitude limits where the r band magnitude is used. We are 
therefore confident that the photometric errors on the SDSS 
LRGs result in an insignificant amount of scattering across 
the selection boundaries. 



4 LUMINOSITY FUNCTIONS 
4.1 Redshift Completeness 

The 2SLAQ survey is not spectroscopically complete unlike 
the SDSS LRG sample which is > 99% complete for input 
target redshifts. Therefore, we must correct the 2SLAQ sur- 
vey for this redshift incompleteness, taking into account any 
dependence on magnitude and/or colour. The redshift com- 
pleteness is defined as the ratio of the number of Sample 8 
2SLAQ LRGs with a reliable redshift to the number of Sam- 
ple 8 target LRGs selected from the SDSS imaging in each 
observed 2dF field. To calculate this, we need to define the 
exact survey area covered by the 2SLAQ survey in order to 
determine the number of possible LRG targets. To achieve 
this, we repeatedly ran the 2dF configuration software on 
random positions until we had configured over 5 million ran- 
dom points in a single 2dF field. This exercise provides a 



detailed map of all possible positions available to the 2dF 
fibres within the field of view. We then built a random cata- 
logue for the whole 2SLAQ survey area by placing this single 
randomised 2dF field at every observed field centre (see Can- 
non et al. 2006). Finally, we remove any regions not in the 
original target input catalogue, i.e., edges of the 2dF fields, 
that extended beyond the SDSS photometry, and holes in 
the SDSS coverage. This produced a random catalogue of 
approximately 400 million positions, covering every possible 
position a 2dF fibre could have been placed throughout the 
whole 2SLAQ survey. We then pixelised these random posi- 
tions (into 30 by 30 arcsecond pixels) to generate a survey 
mask and positively flagged all pixels that contain at least 
one random position. 

The survey mask was used in two ways. First, we calcu- 
lated the area of the survey by summing all positive pixels, 
giving an area of 180.03 deg^. Secondly, we used the mask 
to define those LRGs in the input catalogue that could have 
been included in the 2SLAQ survey in order to calculate the 
redshift completeness. We also restrict the observed LRGs 
in the same manner resulting in about 0.5% of the observed 
LRGs (with redshifts) being excluded. This is caused by 
slight changes made to the 2dF configuration software dur- 
ing the 2SLAQ survey which we are unable to account for 
when constructing our mask. Figure |7| shows the redshift 
completeness of the 2SLAQ survey as a function of mag- 
nitude and colour. There are no significant dependence of 
the redshift completeness on the r — i colour, and i magni- 
tude only shows any significant dependence fainter than i 
= 19.7. However, we do witness a dependency on the g — r 
colour of the LRGs. We correct for this dependence by fit- 
ting a 3rd order polynomial to the data (as shown in Figure 
and use this function to calculate the completeness for 
each LRG depending on its observed g — r colour. Including 
this correction changes the 2SLAQ LRG LF by < 1% com- 
pared to assuming an overall redshift completeness of 76.5% 
regardless of its g — r colour. 
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Figure 5. The colour— magnitude relationship for LRGs satisfying the 2SLAQ selection criteria. We only show here LRGs within 
the 2SLAQ survey area covered by the SDSS multi-epoch photometry. The magnitudes used here were taken from the multi-epoch 
photometry as they are better signal-to-noise (Table 1). The left-hand panel shows LRGs that satisfy the 2SLAQ selection criteria for 
the single— epoch (dashed contours) and multi-epoch photometry (solid contours). The right-hand panel shows LRGs only selected in 
one of the two sets of photometry; red solid contours are LRGs selected in the multi-epoch photometry but missed in the single-epoch 
photometry, and visa-versa for the blue dashed contours. 



Table 2. The Mo. 2^ luminosity function with evolutionary cor- 
rections for the SDSS and 2SLAQ LRGs. The values of Mo. 2^ 
are the median of the bins, the density is given in units of 10~® 



mag ^ 


and the errors are 


1 sigma. 




SDSS 




2SLAQ 


Mo.2r 


Density 


Mo.2r 


Density 


-24.45 


0.03 ± 0.03 


-24.55 


0.20 ± 0.19 


-24.17 


0.15 ± 0.06 


-24.20 


0.32 ± 0.20 


-23.94 


0.62 ± 0.11 


-23.91 


0.63 ± 0.24 


-23.68 


1.91 ± 0.19 


-23.68 


1.99 ± 0.42 


-23.44 


6.51 ± 0.34 


-23.43 


7.72 ± 1.03 


-23.19 


17.56 ± 0.55 


-23.21 


15.45 ± 1.81 


-22.96 


34.47 ± 0.76 


-22.96 


30.59 ± 3.11 


-22.75 


32.59 ± 0.75 


-22.73 


32.25 ± 3.59 


-22.52 


10.47 ± 0.44 


-22.50 


16.94 ± 2.44 


-22.32 


0.86 ± 0.14 


-22.33 


1.19 ± 0.70 



Table 3. The Mo. 2^ luminosity function without evolutionary 
corrections for the SDSS and 2SLAQ LRGs. The values of Mo.2r 
are the median of the bins, the density is given in units of 10~® 
Mpc~^ mag~^ and the errors are 1 sigma. 

SDSS 2SLAQ 



Mo.2r 


Density 


Mo.2r 


Density 






-24.86 


0.20 ± 0.19 






-24.59 


0.14 ± 0.13 


-24.46 


0.03 ± 0.03 


-24.41 


0.45 ± 0.22 


-24.14 


0.23 ± 0.07 


-24.13 


0.78 ± 0.27 


-23.90 


0.76 ± 0.12 


-23.89 


4.50 ± 0.62 


-23.63 


2.84 ± 0.22 


-23.67 


12.18 ± 1.30 


-23.39 


7.96 ± 0.37 


-23.41 


22.29 ± 1.83 


-23.15 


22.15 ± 0.61 


-23.17 


31.95 ± 2.05 


-22.92 


36.64 ± 0.78 


-22.92 


26.79 ± 2.88 


-22.71 


27.68 ± 0.69 


-22.71 


9.54 ± 2.34 


-22.47 


6.96 ± 0.36 







4.2 The Luminosity Function 

In Figure IHl we present the Mo. 2^ luminosity function (LF) 
of the 2SLAQ and SDSS LRG samples as described above, 
both without (top panel) and with (bottom panel) passive 
evolution corrections. We calculated these LFs using the 
1/Vmax method, where for each galaxy we compute the max- 
imum and minimum redshifts at which it would have been 
selected, including the K-l-e corrections described above. For 
all LRGs brighter than Mo.2r = -23, the maximum and min- 
imum redshifts correspond to redshift limits described pre- 
viously, namely 0.17 < z < 0.24 for the SDSS sample and 
0.5 < z < 0.6 for the 2SLAQ sample. Therefore, above this 
absolute magnitude, both samples are essentially volume- 
limited. The luminosity functions were then determined by 
calculating the volume (Vmax) within which each galaxy is 
detectable, modified by the colour-dependent redshift com- 
pleteness corrections for the 2SLAQ sample, and summed 
over all galaxies in the sample. We then apply the photo- 
metric error scattering corrections discussed in Section [3.31 



to each bin. We provide the numerical values of the SDSS 
and 2SLAQ LFs in Tables Hand 01 

We use jack-knife (JK) re-sampling to calculate the 
errors on our luminosity functions. This is achieved by split- 
ting the SDSS and 2SLAQ samples into 20 subregions, of 
equal area, and re-calculating 20 LFs with each of these 
subregions removed in turn. We find that for the 2SLAQ 
sample the JK errors are up to 30% larger than the Poisson 
errors in the four faintest magnitude bins (those which con- 
tain the most galaxies), while the two error estimates are 
the same for the brighter bins. We therefore quote in Tables 
|5|and|3| and plot in Figure |S| the larger of these two errors 
combined with the errors introduced by the photometric er- 
ror scattering correction. For the SDSS sample, we find that 
the JK errors are consistent with the Poisson errors for all 
magnitude bins; we therefore use the Poisson errors. 

Table|l]lists the integrated number and luminosity den- 
sities for a series of magnitude limits for the 2SLAQ and 
SDSS samples and their ratio. These are again corrected for 
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Table 4. The integrated number and luminosity density of the evolution corrected SDSS and 2SLAQ samples and the ratio of these two 
measurements . 



Sample 


Density (xlO"'' Mpc^ 




Luminosity Density (xlO® 


Lq Mpc-3) 




Mo.2r < -22.5 


Mo. 2^ < -23.0 


Mo. 2^ < -23.5 


Mo.2r < -22.5 


Mo. 2,. < -23.0 


Mo. 2^ < -23.5 


SDSS 


25.09 ± 0.33 


9.77 ± 0.20 


1.12 ± 0.07 


3.62 ± 0.05 


1.79 ± 0.04 


0.31 ± 0.02 


2SLAQ 


24.36 ± 1.47 


9.33 ± 0.68 


1.24 ± 0.16 


3.53 ± 0.21 


1.78 ± 0.12 


0.36 ± 0.05 


Ratio 


1.03 ± 0.06 


1.05 ± 0.08 


0.91 ± 0.13 


1.03 ± 0.06 


1.01 ± 0.08 


0.86 ± 0.13 



the effect of the photometric errors scattering galaxies into 
and out of the sample. However, in this instance we use the 
fit to the relation shown in Figure |S| as we need to make a 
correction to each individual galaxy. The errors are again a 
combination of the JK errors and those from the photomet- 
ric scattering. 



5 DISCUSSION 

As seen in Figure|Hl the SDSS and 2SLAQ LFs brighter than 
Mo.2r = -22.6 are in excellent agreement when the passive 
evolution corrections are included. Fainter than this limit we 
are not confident of the photometric scattering correction 
we have made, although we note that the LFs are still in 
reasonable agreement. The agreement of these luminosity 
functions is further confirmed by calculating the integrated 
number and luminosity density of LRGs as given in Table 0| 
Brighter than Mo. 2^ = -22.5, both the integrated luminosity 
and number density of the 2SLAQ and SDSS samples agree 
to within their one sigma errors, and are measured to better 
than 10% out to z = 0.6. 

Throughout the analysis presented herein, we have con- 
sistently used the same simple passive evolution model for 
predicting and correcting the colours and luminosities of 
LRGs as a function of redshift, and this agreement demon- 
strates the lack of any extra evolution, beyond the passive 
fading of old stars, out to z ~ 0.6. This result con firms 
the underlying assum ptions of lEisenstein et all (|20^), and 
ICannon et alJ i2006h . that the majority of LRGs out to 
z ~ 0.6 can be selected via straightforward colour cuts, 
in multicolour data, assuming simple passive evolution of 
th eir stellar popu lations. The result also confirms the work 
of iBernardi et all li2003G ^d'l for lower redshift massive ellip- 
ticals in the SDSS. 

It may appear that our lack of extra evolution beyond 
passive (out to z ~ 0.6) is in confiict wi th recent results 
from the CI 7, DEEP2, and SXI3S sur veys jBell et al.ll20ol 
iFaber et ai]l2005l : lYamada et al.ll2005^ . These smaller-area, 
but deeper (in magnitude limit and redshift), surveys find 
evidence for a change in the density of red galaxies out to 
~ 1 beyond that expected fr om passive fading o f the stel- 
lar populations. For example, iFaber et alJ 1I2OO5I1 report a 
quadrupling of (j)* for red galaxies since « = 1, although 
this result is strongest in their highest redshift bin, where 
they admit their data are weakest. A direct comparison with 
these deeper surveys is difficult because of the differences in 
colour selections used for the surveys, as well as the rela- 
tive luminosity ranges probed by the different surveys, i.e., 
the 2SLAQ survey is designed to probe galaxies brighter 
than a few L* , while the DEEP2, SXDS and C17 surveys 



effectively probe galaxies below L* at z ~ 0.6 (due to their 
smaller areal coverage and fainter magnitude limits). 

However, to facilitate such a comparison, we show in 
FigureEl the LFs from Figure|Hl and the C17 red galaxy LFs 
jBell et al.ll2004 Figure 3) for the same redshift range and 
K-|-e corrected to Mo.2^. We only plot our LFs to Mo. 2^ < 
-22.9 as we do not include all the red galaxies fainter than 
this due to the SDSS LRG selection criteria. Figure|n|demon- 
strates that when one restricts the data to the same red- 
shift range, there is excellent qualitative agreement between 
the 2SLAQ and C17 luminosity functions. We are unable 
to make a quantitative comparison due to the difficulty 
in exactly matching the selection criteria of the two sur- 
veys. Taken together, the surveys shown in Figure |5| ex- 
tend the evidence for no evolution in the LF of LRGs to 
Mo. 2^ < —21, which is close to L* in the LF. Figure El also 
demonstrates that these two surveys are probing different 
luminosity regimes at z < 0.6 as there is at most only 0.5 
magnitudes of overlap in their LFs in which the C17 survey 
is becoming seriously affected by small number statistics due 
to its smaller areal coverage. 

The C17 data presented in Figure 9 agrees with our 
findings that for the brightest galaxies there appears to be 
no evidence for density evolution out to z ~ 0. 6. This re- 
sult is not necessari l y in c o nfiict with the work o f iBell et alJ 
lj2004h : lFaber et alJ |2()0^ 'l : lYamada et alJ lj2005h . as we are 
still probing different redshift and luminosity ranges than 
these other studies. Taken together, these results could in- 
dicate the existence of different evolutionary scenarios above 
and below L* in the luminosity function, i.e., above L* , 
galaxies have only evolved passively (since z = 0.6), while 
below this luminosity, the r ed galaxy population is experi- 
encing significant evolution. iKodama et al.l ^2004) sees sim- 
ilar evidence for differential evolutionary trends with lu- 
minosity at ~ 1, and claim this supports the idea of 
"down-sizing" i.e., the galaxy evolutionary processes (like 
star-formation and assembly) decrease with increased lumi- 
nosity as a function of redshift. These higher redshift obser- 
vations are also consistent with the observed transition, at 
Mr ~ —20.5, in the local colour-magnitude relationship be- 
tween a dominant "red" population of galaxies (above this 
luminosity) , compared t o "blue" population below. Likewise, 
iKauffmann et akl i2003l) find a significant change in the dis- 
tribution of stellar masses of local galaxies at the same lu- 
minosity. A more detailed joint analysis of the 2SLAQ and 
deeper surveys will be presented in a forthcoming paper. 

The luminosity functions given in Figure |H| place tight 
constraints on models of massive galaxy formation and evo- 
lution. Our results appear to favour little, or no, density 
evolution, as we only require the expected passive evolution 
of the luminosities of these LRGs to explain the observed 
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Figure 6. Top and Middle - The magnitude distributions for all 
2SLAQ LRGs in the redshift range 0.5 < 2 < 0.6 with multi- 
epoch photometry. This subsample is further split into LRGs 
selected from the multi-epoch data (dot-dashed line), from the 
single-epoch data (solid line), and those LRGs selected only from 
multi-epoch data but not the single epoch data (dashed line) and 
visa- versa (dotted line). The distributions are normalised to the 
number in the single-epoch subsample. The number of LRGS in 
each of these subsamples is given in parentheses in the top panel. 
The top panel shows those passing the original 2SLAQ selection 
criteria, the middle those passing the additional selection criteria 
discussed in Section 13.21 Bottom - The ratio of number selected 
using the multi-epoch data to the number selected using the sin- 
gle epoch data using the additional selection criteria discussed in 
Section l3.2l The solid line shows a polynomial fit and the dashed 
lines the 1 sigma errors on the fit. 




Figure 7. The redshift completeness of 2SLAQ LRGs as a func- 
tion of apparent i magnitude (top panel), r — i colour (middle 
panel) and g — r colour (bottom panel) . The solid line in the bot- 
tom panel shows the 3rd order polynomial fit to the data and 
used to correct the sample for this incompleteness as a function 
of g — r colour. The bins, chosen to contain 800 LRGs each, are 
plotted at the mean magnitude of the bin with the one sigma 
error bars. 



differences in their LFs as a function of redshift. In other 
words, there are already enough LRGs per unit volume at 
z ~ 0.6 to account for the density of LRGs measured at 
z ~ 0.2. To study this further, we must compare our re- 
sults with the late st predictions f o r mas sive galaxy evolu- 
tion. For example, iDe Lucia et alj ijioO^ have used the ef- 
fects of AGN feedback to regulate new star-formation in 
massive ellipticals within their semi-analytical Cold Dark 
Matter (CD M) model of galaxy form ation. As shown in Fig- 
ures 4 & 5 of iDe Lucia et al.l Jgooi), they find that 50% of 
stars in 2 = massive ellipticals are already formed by a 
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Figure 10. The Mo-2,, luminosity function with passive evolution 
corrections for the 2SLAQ LRGs (solid data points) and fit to 
the SDSS LRGs black line. The lines show the effect of splitting 
varying fractions of the 2SLAQ LRGs in two, simulating major 
mergers between z = 0.6 and z = 0.2. 



FHgure 8. The Mo. 2^ luminosity function without (top panel) 
ami with (bottom panel) passive evolution corrections for both 
the SDSS (open data points) and 2SLAQ (solid data points) LRG 
samples. The points are plotted with their one sigma error bars 
as described in the text. 




Figure 9. The Mo.2r luminosity function with passive evolution 
corrections for the SDSS (blue open data points), 2SLAQ (red 
solid data points) LRG samples, and the COMBO-17 red galaxies 
at z = 0.25 (bla ck open stars) and z = 0.55 (green solid stars) 
jBell et alJl2004) . The dashed lines show the Schechter function 
fit to the COMBO-17 points. The points are plotted with their 
one sigma errors. 



Table 5 . The chi-squared values for fitting difference merger frac- 
tions 



Fraction Mo.2r < -22.75 Mo.2r < -23.25 

Merging Reduced Prob Reduced Prob 



0.68 0.70 0.53 0.76 

0.10 0.64 0.74 0.90 0.47 

0.25 1.66 0.10 2.00 0.07 

0.50 4.54 < 0.0001 3.52 0.003 

0.75 6.07 < 0.0001 6.08 < 0.0001 



median redshift of z — 2.6, yet 50% of the stellar mass of 
z — massive ellipticals is not in place until a median red- 
shift of 2 = 0.8. In Figure 9 of their paper, they show that 
galaxies more massive than ~ 10" M© are built up through 
~ 5 major mergers, that m ust be "dry" (with out gas) to 
prevent new star-formatio n l)van Dokkumll200R|) . Using the 
iBruzual fc Chariot! (l2003l) models described earlier, we es- 
timate that our LRGs have stellar masses > 5 x 10^^ Mq, 
consistent with the massive galaxy sample discussed by 
iDe Lucia et"ai] ^00^. 

We investigate a simple mod el, motivated by t he ide a of 
"dry mergers" and the results of iDe Lucia et all (|200^, to 
simulate the effect on the luminosity function of the hierar- 
chical build-up of these LRGs through major merger events. 
To achieve this, we fit the SDSS LRG LF predict the higher 
redshift (at z — 0.55) 2SLAQ LRG LF under the assumption 
that a given fraction of the SDSS LRGs were formed from 
a major merger of two equal mass progenitors, i.e., we as- 
sume that two 2SLAQ LRGs have merged between z = 0.55 
and z = 0.2 to form a more massive SDSS LRG. We then 
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determine the likely fraction of z — 0.2 LRGs that could 
have been formed this way by fitting (via x^) this model to 
the observed z = 0.55 2SLAQ LRG luminosity function. We 
could have performed this test using the actual data, rather 
than fitting the z = 0.2 SDSS LRG LF, but unfortunately 
such a method would suffer from small number statistics at 
the bright end of the LF, resulting in the brightest bin dis- 
appearing as there are no brighter LRGs being split to refill 
the bin. However, the values are almost identical for the 
other bins. 

The results are presented in Figure [THl and listed in Ta- 
ble [S] They reveal that the 2SLAQ and SDSS LFs are con- 
sistent with each other without any need for merging. At the 
3(7 level, we can exclude merger rates of > 50%, i.e., more 
than half the LRGs at z = 0.2 are already well- assembled, 
with more than half their stellar mass in place, by z ~ 0.6. 
This observation is consistent with Masiedi et al. ( 2003,) who 
find that LRG-LRG mergers can not be responsible for the 
mass growth of LRGs at z < 0.36 based on the small-scale 
clustering amplitude of SDSS LRG correlation function. 

Our limit is barely consistent with the predictions in 
Figure 5 of J)c Lucia et al. ( 20Q5ji) , where they show that 
~ 50% of 2 = massive ellipticals have accreted 50% of their 
stellar mass since z ~ 0.8. We note that our simple model 
does not constra i n the rate of minor mergers; the results of 
iRoseboom et al.) i2006ll on the spectral analysis of 2SLAQ 
LRGs suggests the ~ 1% of our LRGs have experienced a 
small burst of star-formation in the last Gigayear (based on 
the observed H(5 line), which affects less than 10% of their 
stellar mass. 

Our results are a challenge for models of hierarchical 
galaxy formation. More detailed comparisons with semi- 
analytical CDM models are requi red and will be inv esti- 
gated in other papers. For example. iBower et al.l ll2005ll also 
include AGN feedback in their se mi-analytic sim ulations, 
but follow the model suggested bv iBinnevI ll2004ll whereby 
the AGN heating and the gas cooling form a self-regulating 
feedback loop if the gas is in the hydrostatic cooling regime 
(found in groups and cluster) and the central black hol e is 
suitab ly massive. Initial results suggest that the Bower c t alJ 
i2005ri prescription provides a better fit to the LRG evo- 
lution discussed here (Bower, priv. comm.). We also have 
significantly more data than used in this paper, i.e., if we 
could precisely model the K-l-e corrections of these LRGs 
over the joint redshift range of the SDSS & 2SLAQ surveys, 
we would gain a factor of 2 increase in the number of LRGs 
used to compute their luminosity functions. In future work, 
we will investigate t he use of other s tellar synthesis models 
for such corrections 
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APPENDIX A: K AND EVOLUTIONARY 
CORRECTIONS 

We present here tables of the Q — i colours, K and K +e cor- 
rections derived from the lBruzual fc Charlod ll2003ll models 
described in Section ITTI 
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Table Al. K and K+e corrections for the SDSS LRGs to z = 



0.2 assuming the passive model discussed in the text 



K K+e 



z 


9 - « 


u 


9 


r 


t 


u 


9 


r 




0.150 


1.697 


-0.398 


-0.280 


-0.086 


-0.035 


-0.312 


-0.214 


-0.030 


0.017 


0.175 


1.813 


-0.207 


-0.141 


-0.043 


-0.018 


-0.162 


-0.108 


-0.015 


0.008 


0.200 


1.929 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.225 


2.034 


0.202 


0.132 


0.042 


0.017 


0.159 


0.098 


0.015 


-0.008 


0.250 


2.121 


0.422 


0.254 


0.088 


0.043 


0.336 


0.185 


0.032 


-0.008 


0.275 


2.208 


0.647 


0.376 


0.140 


0.069 


0.523 


0.272 


0.056 


-0.007 


0.300 


2.300 


0.879 


0.506 


0.190 


0.097 


0.726 


0.367 


0.079 


-0.005 


0.325 


2.393 


1.116 


0.644 


0.242 


0.128 


0.934 


0.466 


0.104 


0.001 


0.350 


2.450 


1.333 


0.753 


0.303 


0.168 


1.125 


0.537 


0.137 


0.015 



Table A2. K and K-|-e corrections for SDSS LRGs to z = 0.2 assuming the passive plus star-forming model discussed in the text 



K K+e 



z 


9 


— I 


u 


9 


r 




u 


9 


r 




0.150 


1. 


,597 


-0.217 


-0.244 


-0.078 


-0.032 


-0.189 


-0.195 


-0.029 


0.016 


0.175 


1. 


.703 


-0.107 


-0.123 


-0.039 


-0.016 


-0.095 


-0.098 


-0.014 


0.007 


0.200 


1. 


,808 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.225 


1. 


,903 


0.094 


0.113 


0.038 


0.016 


0.086 


0.088 


0.014 


-0.007 


0.250 


1. 


,982 


0.186 


0.217 


0.079 


0.039 


0.174 


0.167 


0.030 


-0.007 


0.275 


2. 


,060 


0.268 


0.321 


0.126 


0.064 


0.261 


0.246 


0.053 


-0.006 


0.300 


2, 


143 


0.340 


0.429 


0.170 


0.089 


0.343 


0.332 


0.074 


-0.003 


0.325 


2, 


,225 


0.402 


0.542 


0.217 


0.118 


0.117 


0.420 


0.097 


0.003 


Q.m) 


2, 


,27,"") 


0. ill 




0.271 


0.1,-),"") 


0. 172 


0.181 


0.12.S 


0.017 



Table A3. K and K-|-e corrections for 2SLAQ LRGs to z = 0.55 assuming the passive model discussed in the text 



K K-l-e 



z 


9-1 


9 


r 


I 


z 


9 


r 


I 


z 


0.450 


2.574 


-0.402 


-0.470 


-0.140 


-0.086 


-0.256 


-0.328 


-0.042 


0.000 


0.475 


2.635 


-0.298 


-0.351 


-0.111 


-0.066 


-0.190 


-0.244 


-0.037 


-0.001 


0.500 


2.682 


-0.204 


-0.227 


-0.075 


-0.049 


-0.131 


-0.157 


-0.025 


-0.006 


0.525 


2.735 


-0.105 


-0.110 


-0.041 


-0.027 


-0.069 


-0.074 


-0.016 


-0.005 


0.550 


2.788 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.575 


2.850 


0.115 


0.103 


0.040 


0.027 


0.078 


0.068 


0.016 


0.005 


0.600 


2.921 


0.245 


0.203 


0.084 


0.058 


0.169 


0.136 


0.036 


0.014 


0.625 


3.001 


0.394 


0.306 


0.137 


0.093 


0.278 


0.208 


0.064 


0.027 


0.650 


3.077 


0.546 


0.414 


0.200 


0.126 


0.389 


0.284 


0.101 


0.037 



Table A4. K and K-|-e corrections for 2SLAQ LRGs to z = 0.55 assuming the passive plus star-forming model discussed in the text. 



K K-(-e 
7, g — i g r i z g r i z 



0.450 


2.376 


-0.244 


-0.412 


-0.126 


-0.080 


-0.189 


-0.308 


-0.039 


-0.001 


0.475 


2.422 


-0.177 


-0.308 


-0.100 


-0.061 


-0.139 


-0.230 


-0.035 


-0.002 


0.500 


2.455 


-0.118 


-0.199 


-0.067 


-0.045 


-0.095 


-0.148 


-0.024 


-0.006 


0.525 


2.491 


-0.060 


-0.096 


-0.037 


-0.025 


-0.049 


-0.070 


-0.015 


-0.005 


0.550 


2.526 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.575 


2.566 


0.063 


0.090 


0.036 


0.025 


0.055 


0.065 


0.015 


0.005 


0.600 


2.609 


0.130 


0.177 


0.076 


0.053 


0.118 


0.129 


0.034 


0.014 


0.625 


2.655 


0.202 


0.268 


0.123 


0.086 


0.191 


0.198 


0.061 


0.026 


0.650 


2.693 


0.269 


0.361 


0.178 


0.117 


0.263 


0.270 


0.095 


0.036 
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Table A5. K and K+e corrections for SDSS LRGs to z = 0.55 



Passive Stajr-forming 
K K+e K K+e 

z g — i u g g r r i u g g r r i g — i u g g r r i u g g r r ^ i 



0.150 


1.697 


0.200 


0.031 


-0.012 





.727 


0.501 


0.386 


1, 


.597 


0.381 


0.067 


-0.004 


0.850 


0.520 


0.387 


0.175 


1.813 


0.391 


0.170 


0.031 


0, 


.877 


0.607 


0.401 


1, 


.703 


0.491 


0.188 


0.035 


0.944 


0.617 


0.402 


0.200 


1.929 


0.598 


0.311 


0.074 


1, 


.039 


0.715 


0.416 


1, 


,808 


0.598 


0.311 


0.074 


1.039 


0.715 


0.416 


0.225 


2.034 


0.800 


0.443 


0.116 


1, 


.198 


0.813 


0.431 


1, 


,903 


0.703 


0.434 


0.118 


1.113 


0.802 


0.429 


0.250 


2.121 


1.020 


0.565 


0.162 


1 


.375 


0.900 


0.448 


1, 


,982 


0.830 


0.573 


0.176 


1.159 


0.876 


0.441 


0.275 


2.208 


1.245 


0.687 


0.214 


1, 


.562 


0.987 


0.472 


2, 


.060 


0.920 


0.685 


0.227 


1.237 


0.954 


0.463 


0.300 


2.300 


1.477 


0.817 


0.264 


1, 


.765 


1.082 


0.495 


2, 


.143 


0.992 


0.793 


0.271 


1.319 


1.040 


0.484 


0.325 


2.393 


1.714 


0.955 


0.316 


1, 


.973 


1.181 


0.520 


2, 


.225 


1.054 


0.906 


0.318 


1.393 


1.128 


0.507 


0.350 


2.450 


1.931 


1.064 


0.377 


2, 


.164 


1.252 


0.553 


2, 


.275 


1.096 


0.995 


0.372 


1.448 


1.192 


0.538 



Table A6. K and K+e corrections for 2SLAQ LRGs to z = 0.2 



Passive Star-forming 
K K+e K K+e 

z g — i r —* g i —* r z —* i r —* g i r z i g — i r ^ g i r z i r —* g i —* r z i 



0.450 


2.574 


-0.781 


-0.214 


-0.183 


-1.043 


-0.458 


-0.408 


2.376 


-0.776 


-0.227 


-0.192 


-1.016 


-0.449 


-0.402 


0.475 


2.635 


-0.662 


-0.185 


-0.163 


-0.959 


-0.453 


-0.409 


2.422 


-0.672 


-0.201 


-0.173 


-0.938 


-0.445 


-0.403 


0.500 


2.682 


-0.538 


-0.149 


-0.146 


-0.872 


-0.441 


-0.414 


2.455 


-0.563 


-0.168 


-0.157 


-0.856 


-0.434 


-0.407 


0.525 


2.735 


-0.421 


-0.115 


-0.124 


-0.789 


-0.432 


-0.413 


2.491 


-0.460 


-0.138 


-0.137 


-0.778 


-0.425 


-0.406 


0.550 


2.788 


-0.311 


-0.074 


-0.097 


-0.715 


-0.416 


-0.408 


2.526 


-0.364 


-0.101 


-0.112 


-0.708 


-0.410 


-0.401 


0.575 


2.850 


-0.208 


-0.034 


-0.070 


-0.647 


-0.400 


-0.403 


2.566 


-0.274 


-0.065 


-0.087 


-0.643 


-0.395 


-0.396 


0.600 


2.921 


-0.108 


0.010 


-0.039 


-0.579 


-0.380 


-0.394 


2.609 


-0.187 


-0.025 


-0.059 


-0.579 


-0.376 


-0.387 


0.625 


3.001 


-0.005 


0.063 


-0.004 


-0.507 


-0.352 


-0.381 


2.655 


-0.096 


0.022 


-0.026 


-0.510 


-0.349 


-0.375 


0.650 


3.077 


0.103 


0.126 


0.029 


-0.431 


-0.315 


-0.371 


2.693 


-0.003 


0.077 


0.005 


-0.438 


-0.315 


-0.365 



